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What is Minimal Model of '^He Adsorbed on Graphite? 
—Importance of Density Fluctuations in 4/7 Registered Solid— 
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Japan 

We show theoretically that the second layer of "^He adsorbed on graphite and solidified 
at 4/7 of the first-layer density is close to the fiuid-solid boundary with substantial density 
fluctuations on the third layer. The solid shows a translational symmetry breaking as in 
charge-ordered insulators of electronic systems. We construct a minimal model beyond the 
multiple-exchange Heisenberg model. An unexpectedly large magnetic field required for the 
measured saturation of magnetization is well explained by the density fiuctuations. The 
emergence of quantum spin liquid is understood from the same mechanism as in the Hubbard 
model and in k-(ET)2Cu2(CN)3 near the Mott transitions. 
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^He layers adsorbed on graphite substrate is a unique two-dimensional correlated Fermion 
system and have continuously offered fundamental issues in condensed matter. In particular, 
adsorption of ^He to the 2nd layer under the corrugation potential of the Ist-layer solid 
shows a variety of phenomena ranging from a correlated Fermi liquid to a solidification at the 
commensurate density of 4/7 relative to the 1st layer. ^ The solid phase bahaves as a quantum 
spin liquid (QSL),-^'^ the nature of which is a long-standing theoretical challenge.^ 

This solidified ^He monolayer has widely been studied by the Heisenberg model with mul- 
tiple spin exchange (MSE).^'^ However, exact diagonalization studies on realistic MSE models 
suggest an opening of spin excitation gap^ in contrast to the gapless nature of QSL revealed 
by specific heat^'^ and magnetic susceptibility measured down to 10 /iK.'^'^ Furthermore, the 
MSE model predicts that the magnetization m saturates above the field /igat ~ 7 Tesla,'''® 
whereas a recent experiment^'^ up to 10 Tesla indicates the saturation at much higher /igat- 

A gapless QSL was reported in numerical studies as the ground state of the two- 
dimensional Hubbard model with geometrical frustration effect near the Mott transition^^"^^ 
supplemented by a report showing numerically the absence of various symmetry breakings. It 
supports the realization of a genuine Mott insulating state without any translational symmetry 
breaking as initially proposed by Anderson.^ This series of studies is indeed relevant and pro- 
vides a realistic model for a subsequently discovered gapless spin liquid in k-(ET)2Cu2(CN)3. 
Although a charge gap exists in Mott insulating states, density fluctuations allowing doubly 
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occupied sites in the Hubbard model near the Mott transition is crucial for the stabilization 
of the QSL. 

However, when we consider the hard core of the interatomic interaction between ^He 
atoms, the Hubbard model with a moderate onsite interaction U near the Mott transition 
with a crucial role of density fluctuations looks unrealistic as a model of ^He monolayer. 

In this letter, we show that the 4/7-density solid is actually located in the vicinity of the 
fluid-solid boundary implying essentially the same character as the QSL found in the Hub- 
bard model with substantial density fluctuations, contrary to the conventional picture. More 
precisely, the density fluctuation in the solid between the 2nd and 3rd layers accompanied by 
a translational symmetry breaking on the 2nd layer solves the puzzles: It causes enhancement 
of the ratio of /isat to the exchange interaction as is revealed in the recent experiment. Fur- 
thermore, it naturally explains why the MSE model is insufficient to describe the 4/7 phase 
of ^He. 

In the 4/7 phase, the 3/4 of ^He atoms on the 2nd layer occupy points just above mid 
points of the edges of triangles formed by the Ist-layer atoms whereas the 1 /4 occupy points 
just above the Ist-layer atoms in a regular fashion as shown in Fig. 1(a). Here, open circles 
represent the atoms on the 1st layer and shaded circles represent actual locations of ^He 
atoms on the 2nd layer when solidified. If '^He atoms are adsorbed on the 1st layer, it forms a 
triangular lattice with the lattice constant a = 3.1826 Aat the saturation density pi = 0.114 
atom/A^.-^ 

The location of the 2nd-laycr atoms is in principle determined as stable points in continuum 
space. In the present treatment, wc simplify the continuum by discretizing it with as much as 
large number of lattice points kept as candidates of the stable points in the solid. To illustrate 
the discretization, we cut out from Fig. 1(a) a parallelogram whose corners are just above 
4 atoms on the 1st layer as in Fig. 1(b). Possible stable locations of '^Hc atoms on the 2nd 
layer are (1) the points just above the mid points of the Ist-layer atoms, (2) the centers of 
the regular triangles and (3) the points just above the Ist-layer atoms. Therefore, we employ 
totally 6 points as the discretized lattice points in a parallelogram as circles in Fig. 1(b). 
Since a unit cell in Fig. 1(a) contains 7 parallelograms, it contains 42 lattice points in total as 
illustrated as circles in Fig. 1(c). Now the 4/7 solid phase is regarded as a regular alignment 
of 4 atoms on 42 available lattice points in the unit cell shown in Fig. 1(c). 

We employ the Lennard-Jones potential 

VLj(r) = 4e[(a/r)^2-(a/r)^], (1) 

for the inter-helium interactionwhere e = 10.2 K and a = 2.56 A.^^ More refined Aziz potential 
is expected to give similar results under this discretization. In the inset of Fig. 2, VLj(r) vs. r 
in the unit of a is shown by the bold solid curve. The interaction term of the lattice model is 
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Fig. 1. (Color online) (a) Lattice structure of the 4/7 phase of •^He. Both Ist-layer atoms (open 
circles) and 2nd-layer atoms (shaded circles) form triangular lattices in the solid phase. The area 
enclosed by the solid line represents the unit cell for the solid of the 2nd-layer atoms (see text). 
The lattice constant of the 1st layer is a. (b) Possible stable location of the 2nd-laycr atoms are 
shown by circles on top of a a x a parallelogram constructed from the 4 neighboring Ist-layer 
atoms, (c) Structure of discretized lattice for the 2nd-layer model. Lattice points are shown by 
circles. 



given by Hy = 'Yliij ('^i is a number operator of a Fermion on the i-th site) with Vij 

taken from the spatial dependence of eq. (1) on the lattice points. In the actual ^He system, 
the chemical potential of the 3rd layer is estimated to be 16 K higher than the 2nd layer. 
^He atoms may fluctuate into the 3rd layer over this chemical potential difference and it is 
signaled by an increase of the specific heat for T > 1 K^' as is reflected by the entropy per 
site larger than k-Q log 2. To take account of this density fluctuation, we here mimic the allowed 
occupation on the 3rd-layer by introducing a simple finite cutoff K:utoff for Vij within the same 
form of Hamiltonian: When ^Ljl^'tj) for fij = |rj — rj| exceeds Kutoff) we take Vij = Kutoff and 
otherwise Vij = l^j(rjj). This allows taking account qualitative but essential part of possible 
occupation on the 3rd layer by the atoms overcoming Fcutoff- We show the case of VcutoS = 16 
K as indicated by an arrow in the inset of Fig. 2. Here, the open circles show V{rij) on the 
lattice sites in Fig. 1(c) for r/a < 2. 

Our Hamiltonian for the lattice model H = /^k + Hy consists of the kinetic energy 
i?K = — Yl{ij){^ij'^\^j + H-C) and Hy. By using the unit-cell index s and the site index / in 
the unit cell, we have rj = + r;. 

After the Fourier transform, q = Cg^i = Xlk Ck.ie''^ ""^/V-^uj the mean-field (MF) approxi- 

-M] 
V 

42 



mation with the diagonal order parameter {n\^^i) leads to H\^H^^ 



l,m=l s' k.p 



(rak,/)np,m - ^(nk,/)(np,m) 



where the inter-atom interaction is expressed as Vij = Vg™ = V {s ) with r^' = rg — rt. Then, 

I 

'V 



we have the MF Hamiltonian Hmf = Hk + H^^. By diagonalizing the 42 x 42 Hamiltonian 
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Fig. 2. (Color online) l^utoff dependence of the "charge gap" Ac. The inset shows the He- He inter- 
action V{r) vs. r (see text). 



matrix for each k, we obtain the energy bands Hmf = X]z=i 

Here we show the results by taking account of the transfers and interactions for \rij\/a < 

2 as indicated by the open circles in Fig. 2. Then, Vij and tij for the ij pairs up to the 
shortest- 19th rjj are retained. For the kinetic energy, several choices of tij are examined and 
here we show the result for tij = to/rfj, assuming that it is proportional to h/{2mrf-). We 
note that the kinetic energy per atom for the 4/7 phase is estimated as 20 K by the path- 
integral Monte Cairo (PIMC) simulation. Hence, we evaluate to by imposing the condition, 
^(ij) *o/rfj{4cj + H.C.)/(4iVu) = 20 K. We thus obtain to = 0.0392 K, which is taken as the 
energy unit. The values of e and a in eq. (1) are given by e/to = 260.14 and a/a = 0.8045, 
respectively. If to is determined so as to reproduce the total kinetic energy of the PIMC result, 
the main result measured in the unit of K shown below is quite insensitive to the choice of 

4. 20 

By solving the MF equations for Hmf, we have the solution of the \/7 x -v/7 commensurate 
structure shown in Fig. 1(a) for Vcutoff ^ 267to = Cutoff - '^^^ "charge gap" opens for Vcutoff > 
Kutoff' ^ shown in Fig. 2. Here, the "charge gap" is defined by Ac = Ef''^^{]i)-E^^^(k), where 
-Ep(k) denotes the minimum or maximum value of the l-th. band from the lowest. The left 
(right) and bottom (top) axes represent Ac and Vcutoff in the unit of K (to), respectively. From 
the specific heat data, Ac is estimated to be ~ 1 K. This leads to Fcutoff/io ~ 300 (namely, 12 
K), which is consistent with the chemical potential difference of the 3rd layer ~ 16 K.^^ Since 
^cutoff /^cutoff ~ ^-^^ ^Z''' phase is located close to the fluid-solid boundary. The effect of 

3 K higher potential on top of the 1st layer ^He than other lattice points of the 2nd layer^^ 
merely shifts the Ac-VcutofT line toward larger Vcutofr: Kutoff i^ changed from ~ 11 K to ~ 14 
K and hence the above conclusion does not change. 

To further understand the nature of the solid near the fluid-solid boundary, we next 
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consider a minimal model H = H}<i + Hu + Hy with ffx = 



H\j = UJ2i n-i^n-ii and Hy = Y.{ij) ^ij J2a.a' '"'ia'Tija', where rij, 



= ct^Cjo- and (ij) denotes the 



pair of the sites. To simulate the quantum phase transition between fluid and commensurate 
solid, we consider a triangular lattice with = 12 sites with A'c = 4 Fermions (see inset 
of Fig. 4). When the nearest neighbor repulsion V = Vij is large in comparison with the 
transfer, a commensurate solid is expected to be realized. To make accurate estimates of 
physical quantities we employ the exact diagonalization. Here the transfer integrals with the 
ath nearest-neighbor ta for a < 3 and the nearest-neighbor repulsion V are retained. We 
take = ^2 = ^3 = 1 and U = V to express the large kinetic energy and the effect of 
Cutoff for '^He. Figure 3 shows the "charge gap". Here, we calculate the ground-state energy 
by introducing the phase factor for the transfer integral: tij = tjjexp[i(^ • (r^ — rj)], where 
(j) = 0ibi -f (j)2^2 with bj being a reciprocal lattice vector which satisfies bj • aj = 5ij. To 
reduce the finite-size effects, the "charge gap" is defined by Ac = max{/x^jj^ — fi^i^x^ 0}, where 
Mmin = min^[^(iVe + 2) - E{N^)]/2 and = nm^4E{N,) - E{N^ - 2)]/2 with E being 

the ground-state energy. We take (j)^ = ^yir/S with ^ = x,y and the integer 7 running from 
to 8, i.e., totally 81 mesh points for iV = 12 and = 18 at the filling n = N^/N = 1/3. 
The results show little system-size dependence indicating the metal-insulator transition at 
= ~ 10 in the bulk limit. The inset in Fig. 3 shows the V dependence of the peak 
value of the equal-time charge and spin correlation functions with the periodic boundary 
condition (b.c), {(pxAy) = (0,0), i.e., iV(q) = I]jjexp[iq- (r^ - r^-)] ((n^nj) - {ni){nj))/N 
andS'(q) = ^-^^^ exp[iq- (r^ -rj)](Sj •Sj)/(3iV). The peak of iV(q) at iqx,qy) = (27r/3, 27r/V3) 
increases rapidly around V = Vc- The peak in S'(q) at {qx, Qy) = (7r/3, 7r/\/3) jumps at a higher 
Vs > Vc suggesting that a commensurate solid for y > is stabilized without a spin order for 
V < Vs implying the QSL for Vc < V < Vg. The realistic choice of V/Vc ~ 1.1 inferred from 
the MF study is located in this QSL region. We note that the ratio of the 3rd-layer promotion 
to the 2nd-layer density is estimated to be about 10-20% from the results of double occupancy 
and the nearest-neighbor {riirij) averaged over the 81 phase factors in = 12 for V/Vc ~ 1.1. 

From the exact diagonalization of = 12 sites with the periodic b.c, the exchange 
interaction J is estimated from high-temperature part of x(^) by the fitting of the high- 
temperature expansion x(^) = (1 — 3J/T)/T on the triangular lattice. By plotting (l/(xr) — 
1)T vs. T as in Fig. 4, we estimate J from the flat part indicated by the arrows. The system- 
size dependence of J is quite small as known in the Hubbard chain, where x(^) high T 
is determined by the local process. Figure 5 (open circle) shows J for each V obtained in this 
way. 

The magnetization is calculated by adding the Zeeman term to H: H — h Sf. We define 
the saturation field /igat at which the total magnetization m = Y^i{Sf) /N reaches its saturation 
value, rrisat = ?^/2. Figure 5 shows /igat for the A'' = 18 sites (open triangle) under the periodic 
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Fig. 3. (Color online) V dependence of "charge gap" on A'^ = 12 (open circle) and N — 18 (filled 
triangle) triangular lattices for ii = ^2 = ^3 = 1 and V — U at n — 1/3. The inset shows V 
dependence of the peak value of iV(q) (filled square) and 5(q) (open triangle) for TV = 12 under 
periodic b.c. 




Fig. 4. (Color online) Temperature and V dependences of susceptibility on = 12 triangular lattice 
at n = 1/3 for ti ^ t2 = ~ 1 and U — V. The inset shows a triangular lattice with = 12 
sites. 



b.c. The present saturation field at V = U = for = 18 reproduces the exact bulk limit 
^sat = 18.0, which is nothing but the width of the n-filled band at /i = 0. This reproduction 
together with slightly smaller hgat for A^ = 12 (see the inset of Fig. 5) suggests that /igat in 
Fig. 5 is close to the bulk-limit. This is one of our central results: /igat and hence hg-^t/J as 
well largely increase in the commensurate solid near the solid- fluid boundary, V = Vc- From 
Fig. 5, we see that 10 Tesla shown as thin lines is still below the saturation magnetic field for 
a realistic choice of V = 11.2 in agreement with the recent experiment. ^'^ 

Although the enhancement of hsat/J also appears at ^2 = ^3 = (not shown), the en- 
hancement is more prominent when t2 is switched on. This is understood by the pertur- 
bation from the large V{= U) limit. When t2 = ^3 = 0, J appears first in the 4th order 
as J(4) = 20tf/(3y^), whereas the 2nd-order term J(2) = ^t^/V and the 3rd-order term 
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Fig. 5. (Color online) V dependence of exchange interaction J (open circle), saturation field /igat 
(open triangle) and its ratio hsat/J (filled square) for ti = t2 = ~ I and V ~ U at n ~ 1/3. 
The inset shows the magnetization process for = 11.2 in iV = 12 (solid bold line) and = 18 
(broken line) under periodic b.c. The difference between iV = 12 and 18 may come from different 
commensurate structures allowed at m ~ nisat/'^- Thin lines (in the inset as well) represent h/J 
corresponding to 10 Tesla in the experiment when we employ J — 0.3 mK.^^ 



J(3) = — 10tft2/y^ appear for t2 ^ 0. For t2 > 0, J(3) becomes ferromagnetic (FM), which 
partially cancels J(4) in J = J(2) + J{3) + ^{4)- This is similar to the cancellation among anti- 
ferromagnetic (AF) J2 > 0, FM — J3 < and AF J4 > in the MSE with n-body exchange 
interactions (— 1)"'J„.^ In short, the enhancement of hsat/J is largely driven by the density 
fluctuations near the fluid-solid boundary, supplemented by the reduction of AF exchange 
through partial cancellation by FM MSE. 

Let us discuss the significance of the density fluctuations near the fluid-solid boundary 
in terms of the observed QSL. The QSL in k-ET2Cu2(CN)3^^ is found in the region of a 
tiny charge gap,^^ which consistently reproduces the QSL numerically found near the metal- 
insulator boundary in the Hubbard model on the triangular lattice.^^'^^'^^ The QSL is sup- 
pressed when the density fluctuations are suppressed at large C/^^ consistently with the absence 
of the QSL phase reported in the spin-1/2 Heisenberg model on the triangular lattice. They 
suggest the importance of density fluctuations for the realization of the QSL in the 4/7 phase 
of ^He. 

The density excitations over the energy Ac ~ 1 K make a peak in the specific heat C(T) at 
T ~ 1 K in addition to a low-temperature peak around T = 10^^ ~ 1 mK reflecting the spin 
excitations, since the exchange interaction is estimated as J(2) = /ycutofr ~ 5x 10^^(tjj7to)^ 
K. The double-peak structure is indeed found in C{T) for y > 10 in the = 12 cluster study 
(not shown) as is observed in ^He.^'^ 

The fluid-solid transition occurs at a very large V^utoff/^o ~ 300 as seen in Fig. 2, which 
reflects the general tendency that the commensurate solid phase dramatically shrinks when 
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the period of the density order becomes long.^^ This explains why the fluid-soUd boundary is 
located near such a large chemical potential difference of the 3rd layer. 

In summary, we have shown that the minimal model for ^He adsorbed on the graphite 
should consider the density fluctuation to the upper layers. In particular, the properties of the 
4/7-solid phase on the 2nd layer are understood only by considering the density fluctuations 
on the 3rd layer, which makes the real system close to the fluid-solid transition beyond the 
description by the MSE model. The magnetic field required for the magnetization saturation 
is largely enhanced in agreement with the experiments. The density fluctuations also serve 
as a key for stabilizing the QSL. Our study predicts that when the lattice constant of the 
Ist-layer solid can be changed, the 4/7 solid phase easily changes to fluid. Experimental tests 
would be highly desired. 
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